Horizontal Gene Transfer as an Indispensable Driver for Evolution of Neocallimastigomycota into a Distinct Gut-Dwelling Fungal Lineage.

Survival and growth of the anaerobic gut fungi (AGF; Neocallimastigomycota) in the herbivorous gut necessitate the possession of multiple abilities absent in other fungal lineages. We hypothesized that horizontal gene transfer (HGT) was instrumental in forging the evolution of AGF into a phylogenetically distinct gut-dwelling fungal lineage. The patterns of HGT were evaluated in the transcriptomes of 27 AGF strains, 22 of which were isolated and sequenced in this study, and 4 AGF genomes broadly covering the breadth of AGF diversity. We identified 277 distinct incidents of HGT in AGF transcriptomes, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. The majority of HGT events were AGF specific (91.7%) and wide (70.8%), indicating their occurrence at early stages of AGF evolution. The acquired genes allowed AGF to expand their substrate utilization range, provided new venues for electron disposal, augmented their biosynthetic capabilities, and facilitated their adaptation to anaerobiosis. The majority of donors were anaerobic fermentative bacteria prevalent in the herbivorous gut. This study strongly indicates that HGT indispensably forged the evolution of AGF as a distinct fungal phylum and provides a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.IMPORTANCE The anaerobic gut fungi (AGF) represent a distinct basal phylum lineage (Neocallimastigomycota) commonly encountered in the rumen and alimentary tracts of herbivores. Survival and growth of anaerobic gut fungi in these anaerobic, eutrophic, and prokaryote-dominated habitats necessitates the acquisition of several traits absent in other fungal lineages. We assess here the role of horizontal gene transfer as a relatively fast mechanism for trait acquisition by the Neocallimastigomycota postsequestration in the herbivorous gut. Analysis of 27 transcriptomes that represent the broad diversity of Neocallimastigomycota identified 277 distinct HGT events, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. These HGT events have allowed AGF to survive in the herbivorous gut by expanding their substrate utilization range, augmenting their biosynthetic pathway, providing new routes for electron disposal by expanding fermentative capacities, and facilitating their adaptation to anaerobiosis. HGT in the AGF is also shown to be mainly a cross-kingdom affair, with the majority of donors belonging to the bacteria. This study represents a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.

Development of an RNA interference (RNAi) gene knockdown protocol in the anaerobic gut fungus Pecoramyces ruminantium strain C1A.

Members of the anaerobic gut fungi (AGF) reside in rumen, hindgut, and feces of ruminant and non-ruminant herbivorous mammals and reptilian herbivores. No protocols for gene insertion, deletion, silencing, or mutation are currently available for the AGF, rendering gene-targeted molecular biological manipulations unfeasible. Here, we developed and optimized an RNA interference (RNAi)-based protocol for targeted gene silencing in the anaerobic gut fungus Pecoramyces ruminantium strain C1A. Analysis of the C1A genome identified genes encoding enzymes required for RNA silencing in fungi (Dicer, Argonaute, Neurospora crassa QDE-3 homolog DNA helicase, Argonaute-interacting protein, and Neurospora crassa QIP homolog exonuclease); and the competency of C1A germinating spores for RNA uptake was confirmed using fluorescently labeled small interfering RNAs (siRNA). Addition of chemically-synthesized siRNAs targeting D-lactate dehydrogenase (ldhD) gene to C1A germinating spores resulted in marked target gene silencing; as evident by significantly lower ldhD transcriptional levels, a marked reduction in the D-LDH specific enzymatic activity in intracellular protein extracts, and a reduction in D-lactate levels accumulating in the culture supernatant. Comparative transcriptomic analysis of untreated versus siRNA-treated cultures identified a few off-target siRNA-mediated gene silencing effects. As well, significant differential up-regulation of the gene encoding NAD-dependent 2-hydroxyacid dehydrogenase (Pfam00389) in siRNA-treated C1A cultures was observed, which could possibly compensate for loss of D-LDH as an electron sink mechanism in C1A. The results demonstrate the feasibility of RNAi in anaerobic fungi, and opens the door for gene silencing-based studies in this fungal clade.

Thermal adaptation strategies of the extremophile bacterium Thermus filiformis based on multi-omics analysis.

Thermus filiformis is an aerobic thermophilic bacterium isolated from a hot spring in New Zealand. The experimental study of the mechanisms of thermal adaptation is important to unveil response strategies of the microorganism to stress. In this study, the main pathways involved on T. filiformis thermoadaptation, as well as, thermozymes with potential biotechnological applications were revealed based on omics approaches. The strategy adopted in this study disclosed that pathways related to the carbohydrate metabolism were affected in response to thermoadaptation. High temperatures triggered oxidative stress, leading to repression of genes involved in glycolysis and the tricarboxylic acid cycle. During heat stress, the glucose metabolism occurred predominantly via the pentose phosphate pathway instead of the glycolysis pathway. Other processes, such as protein degradation, stringent response, and duplication of aminoacyl-tRNA synthetases, were also related to T. filiformis thermoadaptation. The heat-shock response influenced the carotenoid profile of T. filiformis, favoring the synthesis of thermozeaxanthins and thermobiszeaxanthins, which are related to membrane stabilization at high temperatures. Furthermore, antioxidant enzymes correlated with free radical scavenging, including superoxide dismutase, catalase and peroxidase, and metabolites, such as oxaloacetate and α-ketoglutarate, were accumulated at 77 °C.

Isolation and characterization of Rhizobium sp. strain YS-1r that degrades lignin in plant biomass.

AIMS: The aim of this work was to isolate novel lignin-degrading organisms. METHODS AND RESULTS: Several pure cultures of bacteria that degrade lignin were isolated from bacterial consortia developed from decaying biomass. Among the isolates, Rhizobium sp. strain YS-1r (closest relative of Rhizobium petrolearium strain SL-1) was explored for its lignin-degrading ability. Microcosm studies showed that strain YS-1r was able to degrade a variety of lignin monomers, dimers and also native lignin in switchgrass and alfalfa. The isolate demonstrated lignin peroxidase (LiP) activity when grown on alkali lignin, p-anisoin, switchgrass or alfalfa, and only negligible activity was measured in glucose-grown cells suggesting inducible nature of the LiP activity. Analysis of the strain YS-1r genome revealed the presence of a variety of genes that code for various lignin-oxidizing, H2 O2 -producing as well as polysaccharide-hydrolysing enzymes. CONCLUSIONS: This study shows both the genomic and physiological capability of bacteria in the genus Rhizobium to metabolize lignin and lignin-like compounds. This is the first detailed report on the lignocellulose-degrading ability of a Rhizobium species and thus this study expands the role of alpha-proteobacteria in the degradation of lignin. SIGNIFICANCE AND IMPACT OF THE STUDY: The organism's ability to degrade lignin is significant since Rhizobia are widespread in soil, water and plant rhizospheres and some fix atmospheric nitrogen and also have the ability to degrade aromatic hydrocarbons.

Draft genome sequence of Microbacterium oleivorans strain Wellendorf implicates heterotrophic versatility and bioremediation potential.

Microbacterium oleivorans is a predominant member of hydrocarbon-contaminated environments. We here report on the genomic analysis of M. oleivorans strain Wellendorf that was isolated from an indoor door handle. The partial genome of M. oleivorans strain Wellendorf consists of 2,916,870 bp of DNA with 2831 protein-coding genes and 49 RNA genes. The organism appears to be a versatile mesophilic heterotroph potentially capable of hydrolysis a suite of carbohydrates and amino acids. Genomic analysis revealed metabolic versatility with genes involved in the metabolism and transport of glucose, fructose, rhamnose, galactose, xylose, arabinose, alanine, aspartate, asparagine, glutamate, serine, glycine, threonine and cysteine. This is the first detailed analysis of a Microbacterium oleivorans genome.

Draft genome sequence and detailed analysis of Pantoea eucrina strain Russ and implication for opportunistic pathogenesis.

The genus Pantoea is a predominant member of host-associated microbiome. We here report on the genomic analysis of Pantoea eucrina strain Russ that was isolated from a trashcan at Oklahoma State University, Stillwater, OK. The draft genome of Pantoea eucrina strain Russ consists of 3,939,877 bp of DNA with 3704 protein-coding genes and 134 RNA genes. This is the first report of a genome sequence of a member of Pantoea eucrina. Genomic analysis revealed metabolic versatility with genes involved in the metabolism and transport of all amino acids as well as glucose, fructose, mannose, xylose, arabinose and galactose, suggesting the organism is a versatile heterotroph. The genome also encodes an extensive secretory machinery including types I, II, III, IV, and Vb secretion systems, and several genes for pili production including the new usher/chaperone system (pfam 05,229). The implications of these systems for opportunistic pathogenesis are discussed.

Draft genome sequence of Staphylococcus hominis strain Hudgins isolated from human skin implicates metabolic versatility and several virulence determinants.

Staphylococcus hominis is a predominant member of the human skin microbiome. We here report on the genomic analysis of Staphylococcus hominis strain Hudgins that was isolated from the wrist area of human skin. The partial genome assembly of S. hominis Hudgins consists of 2,211,863 bp of DNA with 2174 protein-coding genes and 90 RNA genes. Based on the genomic analysis of KEGG pathways, the organism is expected to be a versatile heterotroph potentially capable of hydrolyzing the sugars glucose, fructose, mannose, and the amino acids alanine, aspartate, glutamate, glycine, threonine, cysteine, methionine, valine, isoleucine, leucine, lysine, arginine, phenylalanine, tyrosine, and tryptophan for energy production through aerobic respiration, with occasional lactate and acetate fermentation. Evidence for poly-gamma glutamate capsule and type IV Com system pili were identified in the genome. Based on COG analysis, the genome of S. hominis Hudgins clusters away from the previously published S. hominis genome ZBW5.

Draft genome sequence of Micrococcus luteus strain O'Kane implicates metabolic versatility and the potential to degrade polyhydroxybutyrates.

Micrococcus luteus is a predominant member of skin microbiome. We here report on the genomic analysis of Micrococcus luteus strain O'Kane that was isolated from an elevator. The partial genome assembly of Micrococcus luteus strain O'Kane is 2.5 Mb with 2256 protein-coding genes and 62 RNA genes. Genomic analysis revealed metabolic versatility with genes involved in the metabolism and transport of glucose, galactose, fructose, mannose, alanine, aspartate, asparagine, glutamate, glutamine, glycine, serine, cysteine, methionine, arginine, proline, histidine, phenylalanine, and fatty acids. Genomic comparison to other M. luteus representatives identified the potential to degrade polyhydroxybutyrates, as well as several antibiotic resistance genes absent from other genomes.

Draft genome sequence of Pseudomonas moraviensis strain Devor implicates metabolic versatility and bioremediation potential.

Pseudomonas moraviensis is a predominant member of soil environments. We here report on the genomic analysis of Pseudomonas moraviensis strain Devor that was isolated from a gate at Oklahoma State University, Stillwater, OK, USA. The partial genome of Pseudomonas moraviensis strain Devor consists of 6016489 bp of DNA with 5290 protein-coding genes and 66 RNA genes. This is the first detailed analysis of a P. moraviensis genome. Genomic analysis revealed metabolic versatility with genes involved in the metabolism and transport of fructose, xylose, mannose and all amino acids with the exception of tryptophan and valine, implying that the organism is a versatile heterotroph. The genome of P. moraviensis strain Devor was rich in transporters and, based on COG analysis, did not cluster closely with P. moraviensis R28-S genome, the only previous report of a P. moraviensis genome with a native mercury resistance plasmid.

Draft Genome Sequences of Five Pseudomonas aeruginosa Clinical Strains Isolated from Sputum Samples from Cystic Fibrosis Patients.

We report here the draft genome sequences of five Pseudomonas aeruginosa isolates obtained from sputum samples from two cystic fibrosis patients with chronic colonization. These closely related strains harbor 225 to 493 genes absent from the P. aeruginosa POA1 genome and contain 178 to 179 virulence factors and 29 to 31 antibiotic resistance genes.

Transcriptomic analysis of lignocellulosic biomass degradation by the anaerobic fungal isolate Orpinomyces sp. strain C1A.

BACKGROUND: Anaerobic fungi reside in the rumen and alimentary tract of herbivores where they play an important role in the digestion of ingested plant biomass. The anaerobic fungal isolate Orpinomyces sp. strain C1A is an efficient biomass degrader, capable of simultaneous saccharification and fermentation of the cellulosic and hemicellulosic fractions in multiple types of lignocellulosic biomass. To understand the mechanistic and regulatory basis of biomass deconstruction in anaerobic fungi, we analyzed the transcriptomic profiles of C1A when grown on four different types of lignocellulosic biomass (alfalfa, energy cane, corn stover, and sorghum) versus a soluble sugar monomer (glucose). RESULTS: A total of 468.2 million reads (70.2 Gb) were generated and assembled into 27,506 distinct transcripts. CAZyme transcripts identified included 385, 246, and 44 transcripts belonging to 44, 13, and 8 different glycoside hydrolases (GH), carbohydrate esterases, and polysaccharide lyases families, respectively. Examination of CAZyme transcriptional patterns indicates that strain C1A constitutively transcribes a high baseline level of CAZyme transcripts on glucose. Although growth on lignocellulosic biomass substrates was associated with a significant increase in transcriptional levels in few GH families, including the highly transcribed GH1 ß-glucosidase, GH6 cellobiohydrolase, and GH9 endoglucanase, the transcriptional levels of the majority of CAZyme families and transcripts were not significantly altered in glucose-grown versus lignocellulosic biomass-grown cultures. Further, strain C1A co-transcribes multiple functionally redundant enzymes for cellulose and hemicellulose saccharification that are mechanistically and structurally distinct. Analysis of fungal dockerin domain-containing transcripts strongly suggests that anaerobic fungal cellulosomes represent distinct catalytic units capable of independently attacking and converting intact plant fibers to sugar monomers. CONCLUSIONS: Collectively, these results demonstrate that strain C1A achieves fast, effective biomass degradation by the simultaneous employment of a wide array of constitutively-transcribed cellulosome-bound and free enzymes with considerable functional overlap. We argue that the utilization of this indiscriminate strategy could be justified by the evolutionary history of anaerobic fungi, as well as their functional role within their natural habitat in the herbivorous gut.

Draft Genome of the Arthrobacter sp. Strain Edens01.

We report the draft genome sequence of Arthrobacter sp. strain Edens01, isolated from a leaf surface of a Rosa hybrid plant as part of the Howard Hughes Medical Institute-funded Student Initiated Microbial Discovery (SIMD) project. The genome has a total size of 3,639,179 bp and contig N50 of 454,897 bp.

The Draft Genome Sequence of Xanthomonas sp. Strain Mitacek01 Expands the Pangenome of a Genus of Plant Pathogens.

We report the draft genome sequence of Xanthomonas sp. strain Mitacek01, isolated from an indoor environment vending machine surface with frequent human use in Stillwater, Oklahoma, USA, as part of the Student-Initiated Microbial Discovery project. The genome has a total size of 3,617,426 bp and a contig N50 of 1,906,967 bp.

Draft Genome Sequence of the Environmental Isolate Chryseobacterium sp. Hurlbut01.

We report here the draft genome sequence of the environmental isolate Chryseobacterium sp. Hurlbut01, isolated from a light switch surface in Stillwater, OK, as part of the Student-Initiated Microbial Discovery (SIMD) project. The genome has a size of 3,899,838 bp and a contig N50 of 321 kb.

Assessing the global phylum level diversity within the bacterial domain: A review.

Microbial ecology is the study of microbes in the natural environment and their interactions with each other. Investigating the nature of microorganisms residing within a specific habitat is an extremely important component of microbial ecology. Such microbial diversity surveys aim to determine the identity, physiological preferences, metabolic capabilities, and genomic features of microbial taxa within a specific ecosystem. A comprehensive review of various aspects of microbial diversity (phylogenetic, functional, and genomic diversities) in the microbial (bacterial, archaeal, and microeukaryotic) world is clearly a daunting task that could not be aptly summarized in a single review. Here, we focus on one aspect of diversity (phylogenetic diversity) in one microbial domain (the Bacteria). We restrict our analysis to the highest taxonomic rank (phylum) and attempt to investigate the extent of global phylum level diversity within the Bacteria. We present a brief historical perspective on the subject and highlight how the adaptation of molecular biological and phylogenetic approaches has greatly expanded our view of global bacterial diversity. We also summarize recent progress toward the discovery of novel bacterial phyla, present evidences that the scope of phylum level diversity in nature has hardly been exhausted, and propose novel approaches that could greatly facilitate the discovery process of novel bacterial phyla within various ecosystems.

Survival of the anaerobic fungus Orpinomyces sp. strain C1A after prolonged air exposure.

Anaerobic fungi are efficient plant biomass degraders and represent promising agents for a variety of biotechnological applications. We evaluated the tolerance of an anaerobic fungal isolate, Orpinomyces sp. strain C1A, to air exposure in liquid media using soluble (cellobiose) and insoluble (dried switchgrass) substrates. Strain C1A grown on cellobiose survived for 11, and 13.5 hours following air exposure when grown under planktonic, and immobilized conditions, respectively. When grown on switchgrass media, strain C1A exhibited significantly enhanced air tolerance and survived for 168 hours. The genome of strain C1A lacked a catalase gene, but contained superoxide dismutase and glutathione peroxidase genes. Real time PCR analysis indicated that superoxide dismutase, but not glutathione peroxidase, exhibits a transient increase in expression level post aeration. Interestingly, the C1A superoxide dismutase gene of strain C1A appears to be most closely related to bacterial SODs, which implies its acquisition from a bacterial donor via cross kingdom horizontal gene transfer during Neocallimastigomycota evolution. We conclude that strain C1A utilizes multiple mechanisms to minimize the deleterious effects of air exposure such as physical protection and the production of oxidative stress enzymes.

The genome of the anaerobic fungus Orpinomyces sp. strain C1A reveals the unique evolutionary history of a remarkable plant biomass degrader.

Anaerobic gut fungi represent a distinct early-branching fungal phylum (Neocallimastigomycota) and reside in the rumen, hindgut, and feces of ruminant and nonruminant herbivores. The genome of an anaerobic fungal isolate, Orpinomyces sp. strain C1A, was sequenced using a combination of Illumina and PacBio single-molecule real-time (SMRT) technologies. The large genome (100.95 Mb, 16,347 genes) displayed extremely low G+C content (17.0%), large noncoding intergenic regions (73.1%), proliferation of microsatellite repeats (4.9%), and multiple gene duplications. Comparative genomic analysis identified multiple genes and pathways that are absent in Dikarya genomes but present in early-branching fungal lineages and/or nonfungal Opisthokonta. These included genes for posttranslational fucosylation, the production of specific intramembrane proteases and extracellular protease inhibitors, the formation of a complete axoneme and intraflagellar trafficking machinery, and a near-complete focal adhesion machinery. Analysis of the lignocellulolytic machinery in the C1A genome revealed an extremely rich repertoire, with evidence of horizontal gene acquisition from multiple bacterial lineages. Experimental analysis indicated that strain C1A is a remarkable biomass degrader, capable of simultaneous saccharification and fermentation of the cellulosic and hemicellulosic fractions in multiple untreated grasses and crop residues examined, with the process significantly enhanced by mild pretreatments. This capability, acquired during its separate evolutionary trajectory in the rumen, along with its resilience and invasiveness compared to prokaryotic anaerobes, renders anaerobic fungi promising agents for consolidated bioprocessing schemes in biofuels production.

Fine-scale bacterial beta diversity within a complex ecosystem (Zodletone Spring, OK, USA): the role of the rare biosphere.

BACKGROUND: The adaptation of pyrosequencing technologies for use in culture-independent diversity surveys allowed for deeper sampling of ecosystems of interest. One extremely well suited area of interest for pyrosequencing-based diversity surveys that has received surprisingly little attention so far, is examining fine scale (e.g. micrometer to millimeter) beta diversity in complex microbial ecosystems. METHODOLOGY/PRINCIPAL FINDINGS: We examined the patterns of fine scale Beta diversity in four adjacent sediment samples (1mm apart) from the source of an anaerobic sulfide and sulfur rich spring (Zodletone spring) in southwestern Oklahoma, USA. Using pyrosequencing, a total of 292,130 16S rRNA gene sequences were obtained. The beta diversity patterns within the four datasets were examined using various qualitative and quantitative similarity indices. Low levels of Beta diversity (high similarity indices) were observed between the four samples at the phylum-level. However, at a putative species (OTU(0.03)) level, higher levels of beta diversity (lower similarity indices) were observed. Further examination of beta diversity patterns within dominant and rare members of the community indicated that at the putative species level, beta diversity is much higher within rare members of the community. Finally, sub-classification of rare members of Zodletone spring community based on patterns of novelty and uniqueness, and further examination of fine scale beta diversity of each of these subgroups indicated that members of the community that are unique, but non novel showed the highest beta diversity within these subgroups of the rare biosphere. CONCLUSIONS/SIGNIFICANCE: The results demonstrate the occurrence of high inter-sample diversity within seemingly identical samples from a complex habitat. We reason that such unexpected diversity should be taken into consideration when exploring gamma diversity of various ecosystems, as well as planning for sequencing-intensive metagenomic surveys of highly complex ecosystems.

Phylogenetic diversity and community structure of anaerobic gut fungi (phylum Neocallimastigomycota) in ruminant and non-ruminant herbivores.

The phylogenetic diversity and community structure of members of the gut anaerobic fungi (AF) (phylum Neocallimastigomycota) were investigated in 30 different herbivore species that belong to 10 different mammalian and reptilian families using the internal transcribed spacer region-1 (ITS-1) ribosomal RNA (rRNA) region as a phylogenetic marker. A total of 267 287 sequences representing all known anaerobic fungal genera were obtained in this study. Sequences affiliated with the genus Piromyces were the most abundant, being encountered in 28 different samples, and representing 36% of the sequences obtained. On the other hand, sequences affiliated with the genera Cyllamyces and Orpinomyces were the least abundant, being encountered in 2, and 8 samples, and representing 0.7%, and 1.1% of the total sequences obtained, respectively. Further, 38.3% of the sequences obtained did not cluster with previously identified genera and formed eight phylogenetically distinct novel anaerobic fungal lineages. Some of these novel lineages were widely distributed (for example NG1 and NG3), whereas others were animal specific, being encountered in only one or two animals (for example NG4, NG6, NG7, and NG8). The impact of various physiological and environmental factors on the diversity and community structure of AF was examined. The results suggest that animal host phylogeny exerts the most significant role on shaping anaerobic fungal diversity and community composition. These results greatly expand the documented global phylogenetic diversity of members of this poorly studied group of fungi that has an important function in initiating plant fiber degradation during fermentative digestion in ruminant and non-ruminant herbivores.